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A model is developed for the confinement and collimation of a baryon poor 
outflow by its surrounding medium. Both, confinement by kinetic pressure of a 
static corona, and confinement by the ram pressure of a supersonic wind ema- 
t— i | nating from a disk surrounding the inner source are considered. Solutions are 

,— 1 1 presented for the structure of the shocked layers of a deflected baryon poor jet 

(BPJ) and exterior wind. The dependence of the opening angle of the BPJ on 
the parameters of the confining medium are carefully examined. It is found that 
the BPJ shock may either converge to the symmetry axis or diverge away from 
it, depending on the opening angle of the BPJ injection cone. In the latter case 
the inner flow exhibits a non-uniform structure, consisting of an ultra-relativistic 
core containing the unshocked BPJ enveloped by the slower, shocked BPJ layer. 
O ' The implications of our results to the prompt GRB emission are briefly discussed. 

o 

Subject headings: gamma rays: bursts - ISM: jets and outflows- MHD - shock 
•£} ' waves 

o : 
r- . 
o ■ 

l> '• 

. !_h 1. Introduction 

Collimation is a generic feature of outflows in astrophysics. It appears in outflows from 
protostars, AGNs, microquasars and GRBs. Despite many efforts to identify the mechanisms 
responsible for the acceleration and collimation of jets, there is as yet no universally accepted 
explanation for these processes. One possibility is that jets are collimated by magnetic hoop 
stress that arises when the field is wound up by rotation. This, however, requires ordered 
magnetic fields of sufficient strength to be present in the vicinity of the central engine. 
Whether such fields can be generated is yet an open issue. Recent 3D simulations of accretion 
disks demonstrate that, although near equipartition magnetic fields can be generated by the 
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MRI process, they are not expected to be ordered in general, except perhaps within a flux 
tube that is connected to the horizon of a rapidly rotating black hole (McKinney & Gammie 
2004; McKinney, 2005; Hawley & Krolik 2006; see also Krolik 2007, and references therein). 
Even if it can be produced, such a field configuration may be kink unstable (Lyubarsky 1992; 
Eichler, 1993; Begelman, 1998; but c.f., McKinney, 2006) and only slowly collimated (e.g., 
Begelman & Li 1994, Lyutikov 2006; but c.f., Beskin & Nokhrina 2006). 

Another possibility is that the relativistic outflow is collimated by the pressure and 
inertia of a surrounding medium. A particular application of the latter mechanism is the 
prompt phase of GRBs, during which a multicomponent outflow is likely to be expelled 
by the central source: The hyper-accretion rates anticipated following the collapse of a 
massive star or coalescence of compact objects render the disk surrounding the black hole 
extremely hot and dense. Under such conditions the disk is optically thick to electromagnetic 
radiation and cools predominantly by emission of MeV neutrinos (Popham et al. 1999; 
Pruet et al. 2003) . The prodigious neutrino luminosity drives a powerful, baryon rich 
wind from the inner disk radii that expands at sub or mildly relativistic speed. The ultra- 
relativistic, baryon poor outflow is produced most likely in a polar region which is protected 
from baryon contamination by either magnetic field lines threading the horizon or by an 
angular momentum barrier. The energy that powers the inner, GRB-producing jet may 
be deposited by neutrino annihilation above the horizon (e.g., Levinson & Eichler 1993; 
MacFadyen & Woosley 1999; Aloy et al. 2000; Rosswog & Ramirez-Ruiz 2003; Aloy et al. 
2005). Alternatively the energy may be extracted magnetically from the spin of the black 
hole (Levinson & Eichler, 1993; Meszaros & Rees, 1997; Van Putten, 2001; McKinney, 2005), 
or the inner disk regions (e.g., Vlahakis & Konigl, 2003; Levinson, 2006a). In this case the 
outflow is expected to be magnetically dominated, however, the magnetic energy should 
eventually be converted somehow to matter and radiation. Collimation of the inner outflow 
by the surrounding wind is expected if the streamlines of the latter diverge slower than the 
streamlines of the former. 

A related issue is the dissipation of the bulk energy of the fireball. The standard 
view (see e.g., Piran 2005, and references therein), at least until recently, was that the 
prompt GRB emission is produced behind internal shocks that form in the coasting region 
of a baryon loaded fireball. This scenario was originally motivated by the detection, in 
many sources, of a power law component of the prompt emission spectrum extending up to 
energies well above the MeV peak, in conflict with the thermal spectrum expected from an 
adiabatically expanding, baryon free fireball (Paczynski, 1986). The post-SWIFT discoveries 
of a shallow afterglow phase at early times, and fastly rising, large amplitude X-ray flares in 
the early afterglow phase (Burrows et al. 2005) introduces new puzzles. According to some 
interpretations (e.g., Fan and Wei 2005) these observations imply prolonged activity of the 
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central engine. If true, then this may indicate that 7-rays are emitted during the prompt 
phase with very high efficiency, in the sense that the remaining bulk energy is a small fraction 
of the total blast wave energy inferred from the afterglow emission. Such episodes can be 
quite naturally explained as resulting from dissipation in an electron-positron dominated 
plasma. Alternative explanations involve variations of the micro-physical parameters and do 
not require extremely high efficiency (Granot et al. 2006). 

The difficulty with purely leptonic fireball models mentioned above may be alleviated 
if dissipation occurs during the outflow acceleration. The claim that internal shocks must 
form in the coasting region, since different shells cannot catch up at radii where the fireball 
is accelerating, applies only to conical geometries. As is shown below, collisions of the 
accelerating fireball with the surrounding wind leads to formation of strong oblique shocks 
well below the coasting region. Bulk Comptonization may then give rise to a nonthermal 
extension of the spectrum well above the thermal peak. It is quite likely that these shocks 
pass through a photosphere in which case the nonthermal photons may escape the system 
before being thermalized. 

A preliminary analysis of hydrodynamic collimation of a baryon poor outflow by a wind 
expelled from a thin torus has been presented in Levinson and Eichler (2000, hereafter LE00). 
For this particular wind geometry LE00 have found that the opening angle equals the ratio 
of the power output of the inner fireball to that of the exterior baryonic wind, and suggested 
that the huge apparent isotropic luminosities implied for some GRB's (e.g., GRB 990123) 
can be understood in terms of collimation of the fireball by an outer baryonic wind with 
more modest energetics, roughly the geometric mean of the fireball energy and its isotropic 
equivalent. Detailed numerical simulations of acceleration and confinement of baryon free 
fireballs produced by deposition of thermal energy above the poles of a non-rotating, stellar 
mass black hole that accretes matter from a thick torus have been performed later by Aloy 
et al. (2005). However, the confining medium in those simulations was restricted to some 
specific configurations obtained in NS+NS and NS+BH merger simulations, which might be 
relevant for short GRBs. Jet breakout in collapsars have been studied using axisymmetric 
hydrodynamic simulations (Aloy et al. 2000). Magnetohydrodynamic simulations of outflow 
formation in collapsars have also been reported recently (Proga, et al. 2003). 

In this paper we generalize the model outlined in LE00 to incorporate additional impor- 
tant features that have been ignored there for the sake of simplicity. In particular, we relax 
the assumption that the baryon poor jet (BPJ) is fully shocked and compute the structure 
of the shocked BPJ layer, explore other geometries of the confining medium, and consider 
also cases where the outer shock has a finite compression ratio. Our approach is to con- 
struct semi-analytic models for the interaction of a baryon poor outflow with a surrounding 
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medium, assuming steady state and axial symmetry. 

2. The model 

We consider an accelerating, ultra-relativistic BPJ confined by the pressure and inertia 
of a surrounding matter. In case of confinement by a supersonic wind, we envision that the 
surrounding wind is expelled (e.g., from an extended disk or torus) over a range of scales 
larger than the characteristic size of the central engine ejecting the BPJ, so that not too 
far out its streamlines may diverge more slowly than the BPJ streamlines, thereby giving 
rise to a collision of the two outflows. In general, such a collision will lead to the formation 
of a contact discontinuity across which the total pressure is continuous, and two oblique 
shocks, one in each fluid, across which the streamlines of the colliding (unshocked) fluids are 
deflected, as shown schematically in fig 1. The details of this structure will depend, quite 
generally, on the parameters of the two outflows and on the boundary conditions. 

Since we are merely interested in the ejection phase following the outflow formation in 
the case of coalescence of a compact binary, or the breakout episode in the case of explosion 
inside a star, as in the collapsar model, it is reasonable to assume that the system under con- 
sideration is time independent and axially symmetric, which greatly simplifies the analysis. 
This of course excludes any temporal effects associated with Kelvin-Helmholtz instabilities, 
mixing of the two fluids, etc., at the interface separating the two outflows which can only 
be studied using 3D numerical simulations. On the other hand, it provides a simple and 
convenient way to study the gross structure of such a system and to gain some insight into 
the physics involved. It may also provide a test case for more sophisticated simulations. To 
this end, we construct in what follows semi-analytic models of colliding outflows. 

2.1. Basic equations 

We seek to determine the structure of the shocked layers assuming that the parame- 
ters of the unshocked fluids are given. The problem is then characterized by the following 
independent variables: the rest mass density pj s , pressure pj s , specific enthalpy Wj S , and 
4-velocity Uj s of the shocked BPJ fluid; the rest mass density p ws , pressure p ws , specific 
enthalpy w ws , and 4-velocity u v ws of the shocked wind fluid, and the cross- sectional radii of 
the inner shock surface rj(z), outer shock surface r w {z) and the contact discontinuity surface 
r c (z). (Henceforth, subscripts j, js, w and ws refer to quantities in the unshocked inner jet, 
in the shocked jet layer, in the surrounding wind, and in the shocked wind layer, respec- 
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tively.) The various zones are indicated in fig 1. The stress-energy tensor associated with 
the inner (outer) shocked layers can be expressed as 

Tjs( ws ) = W js(™s) ul js(ws) ul js(ws) + Pjs(ws)9 flU i (1) 

where is the metric tensor. In the what follows we use cylindrical coordinates (r,<p,z), 
and assume that the system is time independent and axially symmetric. The dynamics of 
the system is governed by mass conservation in each shocked layer: 

d k {p js u%) = 0, (2) 
d k { Pws u k ws ) = 0, (3) 

where the index k in the above equations runs through r and z, energy-momentum conser- 
vation: 

diFjf = o, (4) 

dkT£ = 0, (5) 

and an equation of state for the fluids in each shocked layer. The above system of equations 
is subject to boundary conditions at the shock and contact discontinuity surfaces, obtained 
from the jump conditions across each surface. Continuity of energy, momentum and mass 
across the inner shock surface implies 

C^j's n ju)r=rj — {Tj njy) r=r . } (6) 

{.Pjs^ js^jy) r=r j \Pj^'j^'jv)r=rj > v) 

where nj u = {0,fij) is a space-like unit vector normal to the shock front. Likewise, for the 
outer shock surface we have 



lr=r w ; 



(Pws'U' WS Tl W is^ r=rw (Pw'U' w n W i, S )r=r w - (9) 

Finally, at the contact discontinuity the streamlines of the shocked fluids must be tangent 
to the surface. The jump conditions there then reduce to 

(tf s n cfl ) r=rc = « s n CM ) r=rc = 0, (10) 

(Pjs -Pws)r=r c = 0, (11) 

where n CfJi = (0, n c ) is a unit vector normal to the contact discontinuity surface. 

It is convenient to express the terms on the right hand side of eqs. j6]- [9] explicitly in 
terms of the parameters of the unshocked BPJ and wind fluids, which are given as input, 
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and the impact angles 5j and 5 W (that is, the angle between the velocity of the unshocked 
jet/wind fluid and the corresponding shock surface), which are unknown a priori. Denoting 
by Uj( w ) — ( u j( w ) u j(w)z + u j(«,) u jO)r-) 1//2 the poloidal velocity of the unshocked jet (wind) and 
by lj( w ) — (1 + u2 j{ w )Y^ 2 the corresponding Lorentz factor, the mass and energy fluxes that 
pass through the inner (outer) shock can be written as 

(pj(w)U u j{w) n j{w ) u ) r=rHw) = Pj( w )Uj( w )sm5j( w ), (12) 

( T j(w) n iMv) r=r 3(w) = W j(w )'J j{w) U j{w )Sm5 j{w ). (13) 

The flux of transverse momentum incident into the inner (outer) shock can be obtained by 
projecting the RHS of eq. ([6]) (eq. [8]) on the direction normal to the shock front: 

( T ^) n 3H» n Ku>>)r=r Hw) = w i(w)U 2 j(w) sin 2 5 j{w) + p j(w) . (14) 
Finally, projecting eqs. and ([S]) on the axial direction we obtain 

( T j(w) n j(w>)r=r j(w) = w j[w) u 2 j{w) cos 9 j[w) sin 5 j{w) , (15) 

where 9j( w ) denotes the angle between the velocity of the unshocked BPJ (wind) and the 
symmetry axis at the shock front. The angles 9j and 8 W depend in general on the geometry 
of streamlines. For the BPJ we assume conical streamlines, so that to a good approximation 
tan#j = Tj/z. The impact angle is related to rj(z) through 



(sin , — cos OjV 1 : ) ( r , — /•'■:) 



Sm5j= (l+rf)V2 3 = [(r 2 + ;Z 2 )(1 + r /2 )] l/2' ( 16 ) 

where r'j = drj/dz. Likewise, for the outer wind the relation between 5 W , 6 W and r w is 
obtained upon substituting rj — > r w , 9j — > 9 W and 8j — > — 5 W in eq. ( fT6i) . The angle 9 w (z) is 
computed numerically for each of the flow geometries considered below. 

Under the conditions anticipated the shocked BPJ plasma is radiation dominated. The 
equation of state of the shocked BPJ fluid can be approximated as, 

w js = 4p js + p js c 2 . (17) 

In the shocked wind layer the pressure contributed by baryons may be important. The total 
pressure there is then the sum of the radiation pressure and the pressure contributed by 
protons and electrons: 

Pws = aT* s + 2^kT ws , (18) 
trip 

where T ws is the temperature of the shocked gas, and it is assumed that kT ws /m e c 2 is typically 
less than unity, so that the electrons can be considered non-relativistic. The corresponding 
specific enthalpy is given by, 

w ws = p js c 2 + 4aT^ s + 5—kT ws . (19) 

m p 
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To simplify the analysis further, we suppose that each shocked layer can be approximated 
as a one dimensional, cylindrically symmetric flow along the channel. Specifically, it is 
assumed that inside the shocked layers the flow parameters depend solely on z. Equations 
(J2J)-(j5J) can then be integrated over a section of the flow located between z and z + dz. 
Moreover, as shown in the appendix in the small angle approximation, more precisely to 
order 0(r' 2 ), the poloidal 4- velocity equals the z component of the 4- velocity, viz., Uj s = u z , s 
and likewise for u ws . By employing the latter result we can eliminate two variables thereby 
obtaining a closed set. The details are given in the appendix. Using this procedure, we 
obtain from eqs. (j2j) and (j3J) mass conservation laws for the deflected BPJ and wind fluids: 

d sm o ■ 

-[^-Dhta- f, (20) 

j- [p ws u ws {r 2 w - r 2 c )] = 2p w u w r w Sm — . (21) 

GiZ COS G^uj 

Likewise, energy conservation is obtained from the zeroth component of eqs. (jlj), (j3J): 

— [w Jsl]S u js {r 2 c - r 2 )] = ■IuTjU,!-,^^. (22) 
dz cos ay 

T" [w ws j ws u ws (rl - r 2 )] = 2w w >y w u w r w Sm Sw , (23) 
dz L J cosa^ 

and momentum conservation from the z component of eqs. (jlj) and (J5]): 

d r 9/9 9\ 91 9 dj) j c dv • 

o cos 0i sin 5o . , . 

- [ mi ,u%(rl - r|) - rj ft ,] + rJ-£. + ft^ = 2u> J „?r J — (24) 

C? r 2/2 2\,2 1 2^Vws dr 2 COS ^ sin 5 W / c\ 

— r,,, - r c ) + r„,p ws - r— Pw~r~ = 2w w u w r w . (25) 

dz dz dz cos a™ 

Projecting eqs. (jlj) and (jSj) on the direction normal to the inner and outer shock fronts, 
respectively, we obtain (to order 0(r' 2 ), see appendix) the conditions 

Pj S = WjU 2 sin 2 5j + pj, (26) 
Pws = w w u 2 w sin 2 5 W + p w . (27) 

Finally, eq. (TTTT) implies equal pressures in the two shocked layers, viz., 

Pjs = Pws- (28) 



The system of equations (|T7I) - (l28l) forms a closed set for the 11 independent variables that 
characterize the shocked BPJ and wind fluids. 
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2.2. Input parameters and boundary conditions 

The BPJ is assumed to be injected from a point source into a cone of opening angle 6j )ma x, 
with a uniform energy distribution inside the cone. The exterior wind is envisioned to be 
ejected from a thin disk located at the equatorial plane of the system. We consider separately 
subsonic winds (or coronae), for which confinement is dominated by kinetic pressure, and 
supersonic winds for which kinetic pressure can be ignored. The situation of a static gas 
column with a given pressure profile (e.g., the envelope of a star) is treated as a special case 
of such a subsonic wind. As our canonical model we adopt a split monopole geometry for 
the streamlines of the confining wind (see fig 2), which captures most of the basic features 
in a convenient way, but examine also other wind geometries for comparison. The wind is 
assumed to be uniformly ejected over the range of angles 6 Wy min < &w < tt/2 (in the upper 
half plan), and has a total luminosity L w . The integration starts at the position where 
collision of streamlines first occurs (see fig 2). We find it most convenient to obtain the 
initial values of the flow quantities in the shocked layers using the jump conditions of an 
infinite planar, oblique shock. The direction of the shock surface at the origin is obtained as 
an eigenvalue of the corresponding system of algebraic equations. To check the sensitivity 
of our solutions to the initial values, we made some runs with different initial conditions. In 
most cases it has been found that the solution is highly insensitive to boundary conditions, 
and that after some initial phase the integration quickly converges to the same solution. 



3. Results 

3.1. Pressure confinement 

We consider first the possibility that the BPJ is collimated by the kinetic pressure of the 
surrounding medium. As mentioned above, this is naturally expected when the surrounding 
medium is static or highly subsonic, as in the case where the confining medium is the stellar 
envelope into which a fireball is ejected in the collapsar model for GRBs. This is also expected 
when the confining medium is mildly supersonic if 5 W < -M^ 1 , where M. w is the Mach number 
of the enveloping wind. We envisage that the confining medium has a geometry of a torus. 
To be more concrete, we suppose that within some cylindrical radius R (of the order of a few 
Schwarzchild radii) the region near the injection point of the BPJ is devoid of matter, so that 
collision of the BPJ with the enveloping corona first occurs at a height Zq = Rj tan9j imax 
above the equatorial plan. This choice of boundary conditions merely reflects the expected 
configuration after the initial transient phase (e.g., the corking of the stellar envelope in 
the collapsar model). Formally, the limit of kinetic pressure domination is obtained in our 
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model upon substituting u w = 0, r w = r c in eqs. fl20l) - fl27j) . Equations fl27j) and fl28|) then 
imply p JS = p w . For convenience, we adopt a power law profile for the external pressure, 
p w oc z~ v . When 7] < 4 the oblique shock that forms as a result of collision of the BP J with 
the confining walls should diverge slower than conical or even converge to the axis. This can 
be readily understood by noting that if the shock surface is a cone (i.e., r j = R + tan a jZ 
with atj = const), then the impact angle (eq. [IE] ) decreases as 

sin5 .- (r ? ~ ?L (29) 

at z » i?, and the pressure inside the shocked layer scales as (see eq. [26] ) 

2 2 sin (5 ^' _4 , , 

Pjs ^ • sin J, = — — — - oc z , (30) 

J 2tt(1 - cos6 jtmax )(rj + z 2 ) 

where Lj is the total power and 9j imax is the opening angle of the unshocked jet at the 
injection point. The requirement that the pressure is continuous across the contact discon- 
tinuity surface, Pw = Pj S , implies in turn r\ = 4. Evidently, when 77 < 4 the shock surface 
cannot be conicao An example is shown in fig 3, where results of numerical integration 
are exhibited for p w = p w o(z / zo)~ 3 an d two different values of 0j^ max . In both cases shown 
the contact discontinuity surface has the same shape. The evolution of the inner shock sur- 
face, however, depends on 9j jmax , as can be seen. Whether the shock front converges to the 
axis or diverges away from it depends on the value of the impact angle required to sustain 
momentum balance (see also Komissarov & Falle, 1997). We find that when the condition 
p W Q > Lj sin 2 29j tmax /[16iTcR 2 (l — cos9j jmax )] is satisfied, the shock front is inclined towards 
the axis and the BPJ quickly becomes fully shocked. The reflection of the oblique shocks at 
the symmetry axis should lead to formation of internal shocks, which may affect somewhat 
the resultant structure (e.g., Aloy et al. 2000). In particular it is expected to lead to inho- 
mogeneities in the parameters of the shocked BPJ. However, we believe that this should not 
affect the average structure dramatically. 

In the region where the BPJ can be considered fully shocked we have rj = 0. Equations 
fl22|) and ([24]) then reduce to 

^Pj S lj S Uj S Tir 2 c = Lj, (31) 

and 

^ + l -±U n . = 0, (32) 



1 Using similar arguments Eichler (1982) has shown that if the shocked BPJ fluid cools rapidly, then for 
rj = 4 the cross-sectional radius of the BPJ is independent of z, that is, the BPJ is confined in a cylinder. 
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where Wj s = 4pj S has been adopted. Solving eqs (J5Ij) and fl32|) and using the condition 
Pjs = Pwl- 2 ) we recover the result obtained by LE00, 

r c (z) = r co \p w (z)/p wo \~ 1/4 oc z v/A , (33) 

which holds in the BPJ acceleration zone. Evidently, collimation requires 77 < 4. As pointed 
out in LE00, in the case where the pressure is contributed by a transonic wind, such collima- 
tion will occur naturally in the region located within the acceleration zone of the BPJ and 
above the acceleration zone of the surrounding wind (roughly above its critical point), even 
if the two outflows have the same equation of state. The reason is that in this region the 
density profile of the BPJ declines with radius as r -3 (for a conical geometry) while that of 
the slow wind declines as r -2 , since the wind has reached its terminal velocity. Adopting an 
equation of state of the form p w {n w ) oc for the baryonic fluid yields rj = 2T for a conical 
(or spherical) wind, for which n w (z) oc z~ 2 . Convergence occurs for T < 2, which is the case 
for both relativistic gas (r = 4/3) and non-relativistic gas (r = 5/3). The BPJ breakout 
angle from a collapsar surface can be estimated from equation (1331) to be, 

a c * ~ r' c {z*) = K(r C0 /z )(z*/z ) K ~ 1 , (34) 

where k = rj/4 if the coasting radius lies outside the stellar envelope and k = r]/2 if it lies 
inside, z* ~ 10 10 cm is the location of the stellar envelope (e.g MacFadyen & Woosley, 1999). 
For example in the configuration shown in fig 3 we find that a* = 0.04 rad. 

It is worth noting that collimation shocks should continue to form near the interface 
separating the two outflows whenever deflection of streamlines occurs. Such effects cannot 
be accounted for by the ID model presented here. 

3.2. Confinement by a supersonic wind 

We consider next the situation wherein the BPJ collides with a supersonic wind. The 
analysis can be simplified by assuming that the exterior wind is deflected into a very thin 
shock compared with the cross-sectional radius of the contact discontinuity, and that on all 
scales the momentum transfer is dominated by the ram pressure of the wind. In this limit one 
can take r w = r c in the above equations. This effectively assumes an infinite compression 
ratio for the outer shock and may be a reasonable approximation when the cooling time 
of the shocked fluid is much shorter than the expansion time. However, under conditions 
anticipated in GRBs the optical depth of the shocked wind layer is expected to be very large 
on scales of interest and radiative cooling can be ignored. As is shown below, the resultant 
pressure of the shocked wind layer is larger in this case, owing to the larger inclination of the 
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shock front, thereby giving rise to a better collimation. Nonetheless, for illustrative purposes 
we analyze first the case where the deflected wind layer is negligibly thin. With r c = r w 
equations flTT]), (JUD, fl22D, and (JMD, and the condition 



which follows from eqs. fl26]) - (T28]) . provide a complete description of the BP J structure. 

Further simplification is possible when the inner shock angle is sufficiently small, viz., 
aj << a c , such that the energy and momentum fluxes incident through the inner shock can 
be neglected. Upon taking 6j = 0, the system of equations {T7]), (JUD, (J22D, (ED, and (J32D 
reduce to a single, nonlinear differential equation for the cross-sectional radius r c (z) (LE00): 



Here r c0 = r c {z = zq) and pj S fl = Pj S (z = z ) are the cross-sectional radius and pressure at 
some fiducial point z = z . Solutions of the latter equation can be obtained once the exterior 
wind configuration is specified. LE00 solved eq. (I36p in the case of a wind emanating from 
a thin torus of radius Rt- They have found that at z > Rt the shocked jet is confined 
in a cone of opening angle a c ~ X~ l i where x = L w /Lj is the ratio of the wind and jet 
luminosities. For the split monopole wind depicted in fig 2 we have tan 6^ = r c /(z + d) and 
w wU^Atcc cos 8 ' Wt min[ r 1 + ( z + d) 2 ] = L w . Equation f[3"6"j) then yields r c = r c0 + tana c0 (2 — z ) 
at z > Zq » d, where a c0 denotes the value of a c (z) at z = z , with some nontrivial 
dependence of a c o on x- Of course a c o is limited by geometry, that is, a c o cannot be smaller 
than 6 w ^min- It can be readily shown that this behaviour is not limited to the two examples 
discussed above; for any wind configuration that satisfies t&n9 w ~ r c /z and WwU^ oc z~ 2 
at large z, the shocked jet becomes conical, with an opening angle that depends on the 
geometry and power of the unshocked wind. 

The above trend have been confirmed by numerical integration of the full set of equations 
f|T7|) . (12"U1) . f[2"2l . fliM|) . and . An example is shown in fig 4 for the split monopole wind 
exhibited in fig 2 with 9 w ^ m i n = 8.5°. The inner jet in this example has an initial opening 
angle of 6j >m ax = 45°, and its power is normalized such that Lj/4ircR 2 = 2 x 10 26 ergs cm -3 , 
where R = d tan Q w ^ m in is the radius at which the innermost streamline of the split monopole 
wind intersects the equatorial plane (see fig 2). Collimation occurs typically when x exceeds 
unity. We find that for values of x n °t too large (x < 200 for the choice of parameters in 
fig 4) the inner shock diverges away from the axis. The coaxial structure that forms then 
consists of an inner, fast core (the unshocked jet), and a slower sheath surrounding it (the 
shocked layer), as seen in the left panels in fig 4. For very large values of % the shock surface 
converges towards the axis and the inner jet becomes fully shocked, as seen in the right 



WjU 2 sin 2 Sj + Pj = w w u 2 w sin 2 5 W + p, 



(35) 



(r c /r c0 ) 



-l 



w w u 2 w (sinfl w - cos 6 w r' c ) 2 
Pjs,o (1 + r' 2 ) 



(36) 
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panels in fig 4. At sufficiently large radii the contact discontinuity surface indeed becomes 
conical as expected, with an opening angle that scales roughly as a c0 oc x 1 ^ 3 a ^ « c o > 9 w ,min- 
The dependence of a c0 on the power ratio x is presented in fig 5. The above results depend 
only weakly on the opening angle of the BP J injection cone 9j^ max . Changing the latter from 
45° to 60° altered the above results by no more than a few percent. 

A stronger dependence of a C Q on x has been found for exterior wind configurations in 
which some streamlines are inclined towards the axis, as in the case of a wind from a thin 
torus invoked in LE00, where a c0 oc This stronger dependence is a consequence of the 
large impact angles 5 W that give rise to a much stronger compression near the injection point 
of the BPJ. 

Better collimation can be achieved when the streamlines of the exterior wind diverge 
slower than conical. An example is a parabolic wind, for which the family of streamlines is 
defined by the equation 

\l/(r, z) = V z 2 + r 2 — z = r , (37) 

with ro being the family parameter. A particular choice of ro corresponds to the radius at 
which a given streamline meets the disk. The impact angle is given by 

Tq \/ Z 2 -j- T 2 — Z 

tan^ = — = , (38) 

r r 

which at z » r satisfies, tan 6^ — > r/(2z). The power enclosed between two neighboring 
surfaces \1/ = r and ^ = r + cIvq: 

2 2 2 V drn + 2dr z 
dL = w w u w 2irrdr cos 9 W = w w u w 2nr — (39) 

is conserved. At z » r the last equation yields dL ~ w w v? w inz \fr$ dr , implying u^u^ oc 
z~ l . Assuming again that the BPJ is fully shocked, substituting the above results into eq. 
( |36l) and solving for r c (z), one readily finds r c oc 7 oc z 1 ^ 2 , and pj S oc z~ 2 . Integration of 
the full set ( 1T71) . ( 1201) . ( 1221) . f l24"l) and ( 135]) confirms again this behaviour. Because of the 
larger flux of transverse momentum through the outer shock, convergence of the inner shock 
towards the axis occurs at much lower values of x than in the split monopole case. 

We now consider the limit where energy losses by the deflected wind layer are negligible. 
In that case the deflected wind layer cannot be considered thin and the constraint r w = r c 
must be relaxed. We must then numerically integrate the full set of equations (|T7|) - (|28|) . 
As explained above, eqs. (|T7j) - (l28l) have been derived using the small angle approximation. 
However, for the cases examined below we find that the angle of the outer shock surface, a w , 
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is typically not small. This leads to an overestimate of the ratio of energy and momentum 
fluxes incident through the shock front, and tends to produce some artificial dissipation inside 
the shocked wind layer. To correct for this effect, we computed the energy and momentum 
fluxes that pass through the shock locally by solving the jump conditions of an infinite 
planar oblique shock. Results of such calculations, using for a comparison a split monopole 
wind with the same parameters used in fig 4, are exhibited in fig 6. Comparison of figs 4 
and 6 reveals that better collimation is achieved in the absence of energy losses behind the 
shock. The BPJ in this case does not become conical asymptotically, as in the thin shock 
case exhibited in fig 6, but rather the profile of the contact discontinuity surface approaches 
r c oc z 2//3 roughly. The reason for this behaviour is that the smaller compression ratio of the 
shock gives rise to a larger inclination of the shock front and, hence, a larger impact angle 
S w , as seen in fig 6. This in turn leads to a larger pressure behind the shock, as implied by 
eq. (EZD. 

Our analysis does not take into account effects associated with the formation of a rar- 
efaction wave near the interface separating the shocked BPJ and wind fluids in the region 
where the deflected wind layer is highly supersonic. This should give rise to a somewhat 
steeper decline of the pressure supporting the BPJ, and the consequent alteration of the 
resultant structure in this region. Under certain conditions the fluid that resides between 
the contact discontinuity and the rarefaction wave can quickly accelerate to high Lorentz 
factors (Rezzolla et al. 2003; Aloy & Rezzolla 2006). We find that the contact discontinuity 
surface becomes causally disconnected from the outer shock front ( that is, the shocked wind 
layer expands sideways at a speed that exceeds the local sound speed) at z/R > a few. At 
this height the BPJ has already been substantially collimated. For x = 18 and x = 576 we 
find a c ~ 12° and 6°, respectively at the decoupling point. 



4. Limitations of the Model 

The main limitation of our model is that it cannot account for any temporal effects, par- 
ticularly those associated with Kelvin-Helmholtz instabilities at the interface separating the 
two outflows. Linear analysis by Hardee & Hughes (2003) for a cylindrical, non accelerating 
flow embedded in an external medium indicates that non-axisymmetric modes, particularly 
low order helical and elliptical modes, are likely do grow on time scales of the order of the 
jet crossing time. The presence of a viscous boundary layer (that they modeled as a wind) 
leads to suppression of the non-axisymmetric modes and an increase of the growth rate of 
the zeroth order axisymmetric pinch mode. If the instability grows to a nonlinear state over 
the expansion time then pinching of the fast jet and mixing of wind material with the light 
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jet fluid near the interface is anticipated. This can modify the transverse structure of the 
jet, and produce internal shocks that can lead to further dissipation of the bulk energy of the 
jet (see next section for further discussion). Numerical simulations exhibit some evidence 
for such instabilities (e.g., Aloy, 2005; Hardee & Hughes 2003), although the details may 
depend on the configuration of the confining medium, and on the structure of the viscous 
boundary layer. In particular, we are not aware of any simulations of jet confinement by an 
external, supersonic wind. We naively expect suppression of the instability in cases where 
the external wind is mildly relativistic, owing to the small velocity shear. 



5. Summary & conclusions 

We have constructed a class of semi-analytic models for the collimation and confinement 
of an accelerating, ultra-relativistic outflow by the pressure and inertia of a surrounding 
medium, under conditions anticipated near the central engines of GRBs. We considered two 
different scenarios: In the first one the inner flow is assumed to be confined by the kinetic 
pressure of a subsonic (or static) corona that extends vertically over several characteristic 
scales. In the second one it is envisaged that the inner outflow collides with a supersonic, 
baryon rich wind emanating from a disk or a torus surrounding the inner source. Confine- 
ment in this case is accomplished by the ram pressure of the exterior wind. In general, the 
collision of the inner jet with the external medium creates a strong oblique shock across 
which the streamlines of the jet are deflected. In the case where the inner jet is confined 
by the ram pressure of a supersonic wind a second shock forms across which the exterior 
wind gets deflected. The shocked jet and wind layers are separated by a contact disconti- 
nuity surface across which the pressure is continuous. The gross features described above 
are observed in numerical simulations (Alloy et al. 2000, 2005), although temporal effects 
complicate the structure. The evolution of the shocked jet layer depends, quite generally, on 
the configuration of the confining medium (the pressure profile in the case of a static corona 
or the wind geometry in the case of confinement by a supersonic wind), and the opening 
angle of the injection cone of the inner outflow. We identified cases where the inner shock 
converges towards the axis, and cases where it diverges away from the axis (depending on 
the opening angle of the inner flow). In the former case the inner flow quickly becomes fully 
shocked, whereas in the latter case the structure of the polar outflow consists of a core con- 
taining the unshocked ultrarelativistic jet enveloped by the shocked jet layer that expands 
relativistically, but with a Lorentz factor considerably smaller than that of the unshocked jet. 
The observational consequences of such a structure have been discussed elsewhere (Eichler 
and Levinson 2003, 2004; Levinson and Eichler 2004, 2005). 
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The structure of the inner jet, when confined by the ram pressure of a supersonic wind, 
depends on the geometry of the wind and the conditions inside the shocked wind layer. For 
wind geometries with radially diverging streamlines (e.g., split monopole) we find that if the 
shocked wind layer can be considered negligibly thin (e.g., due to rapid cooling), then after 
an initial collimation phase the inner jet becomes conical with an opening angle that scales 
roughly as (Lj/L^) 1 / 3 above the geometrical limit, where Lj and L w are the total power of 
the unshocked jet and wind, respectively. A stronger dependence of the opening angle on the 
power output ratio has been found for exterior wind configurations in which some streamlines 
are inclined towards the axis, as in the case of a wind from a thin torus invoked in LEOO. 
In the absence of energy losses from the shocked wind layer the streamlines of the shocked 
baryon poor flow diverge slower than conical, owing to the larger pressure of the shocked 
wind. For the cases examined above we find that the cross sectional radius of the contact 
discontinuity scales roughly as r c oc z 2 ^ 3 . However, at some radius the outer shock and 
the contact discontinuity become causally disconnected, at which point collimation should 
be less effective. The opening angle of the inner jet at the decoupling point is typically 
smaller than that found in the case of an infinite compression ratio shock. The formation 
of a rarefaction wave may act as a booster on the BPJ fluid contained between the contact 
discontinuity and the rarefaction wave, as demonstrated by Aloy & Rezzolla 2006, although 
this effect should be explored further for the situation considered here. Free neutrons leaking 
across the interface separating the baryonic wind and the shocked BPJ layer will be picked 
up by the BPJ, leading to baryon loading of the inner flow and the consequent emission 
of high-energy neutrinos and 7-rays, as discussed in Levinson & Eichler (2003). However, 
Kelvin- Helmholtz instabilities at the interface may alter this picture in ways that need to be 
explored. 

In situations in which the inner jet is confined by the kinetic pressure of a medium having 
a pressure profile of the form p oc z~ v , the cross-sectional radius of the shocked baryon poor 
jet scales as r c oc z 71 ' 4 . If the pressure is contributed by a coasting wind, then collimation 
(77 < 4) is expected to occur naturally in the acceleration zone of the inner, baryon poor 
flow. 

It has been argued that internal shocks cannot form in the acceleration zone of a fireball 
since different fluid shells cannot catch up there. As demonstrated above, collisions of the 
accelerating baryon poor outflow with exterior medium can lead to dissipation of a consider- 
able fraction of its bulk energy below the coasting region through formation of strong oblique 
shocks. Moreover, reflection of those shocks at the symmetry axis, which may occur under 
certain conditions, and/or intermittent ejection of the inner outflow and/or the confining 
wind should lead to large temporal variations of the Lorentz factor of the deflected flow at 
a given position, and the consequent collisions of different blobs of baryon poor fluid that 
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have been deflected by the confining medium (see also Eichler 1994; Aloy et al. 2000). The 
internal shocks thereby created would give rise to additional dissipation over a large range 
of radii. Those shocks are expected to be radiation mediated. The radiative efficiency is 
expected to be very high if the photosphere is located below the coasting region. The photon 
distribution produced inside those shocks is likely to have a nonthermal extension at high 
energies owing to multiple Compton scattering on the converging flow (a relativistic ver- 
sion of the Blandford Pyne process). If the shocks are produced at a modest optical depth 
then the nonthermal photons may escape before being thermalized, giving rise to a non- 
thermal spectral component of the prompt GRB emission, as proposed earlier (e.g., Eichler 
1994; Levinson 2006b, and references therein). Detailed calculations are required to assess 
whether this process can indeed explain the observed spectra. 

Finally, we comment on the applicability of this model to other systems. Opacity 
arguments suggest that the rapid variability of the TeV emission observed in the TeV blazars 
implies high Doppler factors of the emitting plasma in the 7-ray emission zone (e.g., Levinson 
2006b). Such high values of the Doppler factor are in clear disagreement with the much 
lower values inferred from radio observations (e.g., Urry & Padovani 1991; Marscher 1999; 
Jorstad et al. 2001). Various explanations, including jet deceleration (Georganopoulos 
& Kazanas 2003), a structure consisting of interacting spine and sheath (Ghisellini et al. 
2005), and opening angle effects (Gopal-Krishna et al. 2004) have been proposed in order to 
resolve this discrepancy. Stationary radio features have been observed also in other systems, 
and on various scales, where relativistic expansion is inferred or expected (e.g., HST1 knot 
in M87). We emphasize that, while the Lorentz factors inferred from radio observations 
reflect the speed of the emission pattern, the rapid variability of the very high energy 7-ray 
emission constrains the Lorentz factor of the fluid elements emitting the 7 rays. Dissipation 
behind recollimation shocks can naturally account for large differences between the speeds 
of the dissipation pattern and the outflowing matter. In particular, stationary features in a 
relativistic outflow can be easily produced, as demonstrated by Komissarov & Falle (1997). 
While it is quite likely that the TeV flares in the BL Lac objects are produced by decelerating 
fronts, recollimation on small scales, particularly in the case of efficient radiative cooling, 
can alleviate the requirements on the Lorentz factors (Levinson 2007, in preparation). It 
has been proposed recently (Stawarz et al. 2007) that the HST1 knot in the M87 jet may 
be associated with a converging shock in a reconfinement nozzle. 

We thank M. Aloy for a careful reading of the manuscript and for constructive criticism. 
This work was supported by an ISF grant for the Israeli Center for High Energy Astrophysics. 
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A. Derivation of the flow equations 

We divide the shocked layers into slices of thickness dz along the symmetry axis (see fig 
7) . Equations (j2J) and (jlj) are then integrated over a volume of the shocked BP J enclosed by 
the inner shock surface, the contact discontinuity surface, and two planes perpendicular to 
the z-axis, located at the points z and z + dz. Using Gauss's theorem we obtain, 
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where hj and h c are the normals to the inner shock and contact discontinuity surfaces, respec- 
tively, dflj and dfl c denote integration over a portion of the shock and contact discontinuity 
surfaces enclosed between z and z + dz, respectively, and dSj and dS c are the corresponding 
surface area elements, as indicated in fig 7. The latter are given explicitly by 
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where the angles a 3 - and a c are related to the cross- sectional radii Tj(z) and r c (z) through: 
drji c \/dz. Recalling that the BP J streamlines are tangent to the contact dis- 
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continuity surface we have u^ s h ck 



T^n ck 



and Tgh ck 



Pj S sin a c . Now, the energy 
flux TV* = Wj S r yj S Uj S involves the poloidal velocity Uj S that appears in the Lorentz factor 
7 JS = (u 2 j S + l) 1 ^ 2 , which on any given streamline of the shocked BPJ fluid is given in terms 
of the angle between the direction of the streamline and the axis, if)(r, z), as: Uj s = Uj S cosi/j. 
Since ip < a c we have 2(1 — Uj S /uj s ) < a 2 c < (dr c /dz) 2 at small angles. To order (dr c /dz) 2 
we can therefore replace it* by Uj thereby eliminating one variable. Under the assumption 
of our ID model, that the flow quantities depend only on the coordinate z, the above results 
yield, 
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By employing the jump conditions at the shock front and the contact discontinuity surface 
(eqs. [6] and [7j), and using eqs. f|T2l and ([TBI we finally arrive at eqs fl20|) . fl22|) and fl24|) . 
Likewise, eqs. ()3]) and (jSJ) can be integrated over the volume of the shocked wind enclosed 
between z and z + dz. Using the same procedure as above we obtain eqs. (|2T|) . (j23"j) and fl25l) 
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Two additional relations can be obtained by projecting the jump conditions (jSJ) and (jSJ) 
on the directions normal to the inner and outer shock surfaces, respectively, using eqs. (pQ) 
and flUD: 



The terms associated with the ram pressure of the downstream fluid (the second terms on 
the left-hand side of eqs. |A6j and |A7j ) depend on the direction of deflected streamlines 
just behind the shock surfaces, which cannot be computed under the assumptions of our 
model. However, these terms are typically smaller than the post shock kinetic pressure. For 
a planar shock w s (u k s fik) 2 jp s = 1/6 for a non-relativistic, radiation dominated shock, and 
w s {u k s hk) 2 /ps = 1/2 in the extreme relativistic limit. For an oblique shock this ratio is even 
smaller. Thus, those terms can be neglected, whereby eqs. (1261) and ( 1271) are obtained. 
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Fig. 1. — A sketch of the colliding outflows model. The unshocked BPJ, shocked BPJ, unshocked 
wind, and shocked wind zones are designated as j, js, w, and ws, respectively. 
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Fig. 2. — A sketch of the split monopole wind geometry. The streamlines of the unshocked wind 
are represented by the dot-dashed lines. The m illinium angle of the wind. W rnin 

is also indicated. 

The dashed and solid lines mark the streamlines of the unshocked BPJ and the inner and outer 
shock surfaces, respectively. The BPJ is assumed to be injected into a cone of opening angle Oj, ma x, 
as indicated. 
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Fig. 3. — Lorentz factor (top) and logarithm of the temperature (bottom) of the shocked and 
unshocked BP J in the case of confinement by a static corona with a pressure profile p w oc z~ 3 , for 
two different angles of the BP J injection cone, as indicated. The dashed line marks the inner shock 
surface. 
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Fig. 4. — Lorentz factor (top) and logarithm of the temperature (bottom) of the shocked and 
unshocked BP J in the case of confinement by a split monopole wind with 6 Wt mi n = 8.5°, for two 
values of the luminosity ratio x, as indicated. In both cases shown Lj/AircR 2 = 2 x 10 26 erg cm -3 , 
where R = d tan 9 W) min is the radius at which the innermost streamline of the wind intersects 
the equatorial plane. The blue solid lines represent the streamlines of the unshocked wind. The 
shocked wind layer is assumed to be negligibly thin in this calculation (that is, the outer shock front 
coincides with the contact discontinuity surface). The opening angle of the BP J injection cone in 
both cases shown is 9j^ max = 45°. The opening angle of the shocked BP J (specifically, the contact 
discontinuity surface) approaches 27° and 13° at z/R ~ a few for x = 18 and \ = 540, respectively. 
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Fig. 5. — Dependence of a C Q on x, computed for a split monopole wind with Q w , m in = 8.5°. 
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Fig. 6. — Same as fig 4, but with no energy losses from the shocked wind layer. The parameters of 
the split monopole wind and the angle 0j tTnax are the same as in fig 4. The inner shock is marked 
by the dashed line and the contact discontinuity surface by the thick solid line. The colored region 
to the right of the thick solid line corresponds to the shocked wind layer. 
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Fig. 7. — Schematic diagram showing the slice of shocked BPJ and wind layers, and the cor- 
responding contours of integration. The various quantities used in the derivation of the model 
equations are indicated 



